Controlling the density of the 2DEG at the SrTiOs/LaAlOs interface 
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The polar discontinuity at the SrTiOa/LaAlOs interlace (STO/LAO) can in principle sustain an 
electron density ol S.SxlO'^* cm~^ (0.5 electrons per unit cell). However, experimentally observed 
densities are more than an order of magnitude lower. Using a combination of first-principles and 
Schrodinger-Poisson simulations we show that the problem lies in the asymmetric nature of the 
structure, i.e., the inability to form a second LAO/STO interface that is a mirror image of the first, 
or to fully passivate the LAO surface. Our insights apply to oxide interfaces in general, explaining 
for instance why the SrTiOs/GdTiOs interface has been found to exhibit the full density of 3.3 x 10^* 
cm~^. 

PACS numbers: 73.20.-r,73.40.Lq,73.61.Le 



The realization of a two-dimensional electron gas 
(2DEG) at the SrTiOs/LaAlOa interface (STO/LAO) 
has set off an explosion of interest in oxide electronics. 
This 2DEG exhibits densities that are difficult to achieve 
in conventional semiconductor&ir— and displays unique 
behavior including ferromagnetism,* superconductivity^ 
and even the puzzhng coexistence of both.~ It has been 
proposed as the basis for novel electronic devices that 
exploit strong electron-electron correlation in the nar- 
row bands derived from d states of the transition metal. '' 
STO/LAO heterostructures have been fabricated using 
methods that allow unprecedented control over layer 
thickness, such as pulsed layer deposition and molecu- 
lar beam epitaxy. ^'■^ While great progress has been made 
in characterization and exploitation of the physical phe- 
nomena, the mechanisms that determine the density of 
electrons in the 2DEG have remained a subject of intense 
debate^"— This lack of understanding inhibits achieving 
the control that is required for device applications. 

Typically, an LAO layer of less than 20 nm is deposited 
on a Ti02-terminated [001]-oriented STO substrate or 
epilayerj^iS, as shown in Fig. [T] The electrical conductiv- 
ity at the buried interface is then probed as a function of 
temperature i^i^ Carrier densities up to 2 x 10^'^ cm~^ have 
been reported.^^^"— Low-temperature sheet resistance 
varying from 10~^ to 10"* fi/D has been observed, dis- 
playing a strong dependence on the oxygen partial pres- 
sure in the growth environment or post-growth annealing 
treatments A dependence of the sheet resistivity 
and carrier density on the thickness of the LAO layer 
has also been observed: heterostructures with LAO lay- 
ers less than 4 unit cells ('^1.6 nm) thick exhibit insulat- 
ing behavior, while thicker layers become conductingfiZ 
However, experiments by Huijben et al}'^ indicated the 
existence of conducting interfaces for LAO layers as thin 
as two unit cells in the presence of an STO layer on top 
of the LAO. 

A variety of models have been put forth to explain the 
origin of the carriers at the interface. ^^^^ However, none 
can account for all the experimental observations, raising 
questions about their validity and general applicability. 
The prevailing models assume that the carriers originate 
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FIG. 1: (Color online) Layer structure of an STO/LAO het- 
erostructure with Ti02-LaO planes at the interface. Nominal 
charges are indicated above each layer. LaO planes act as 
electron donors; the Ti02 plane terminating STO is charge 
neutral, and therefore the interfacial LaO plane acts as a 
delta-doped layer of donors with a density of 0.5 electrons 
per unit cell. 



from somewhere other than the ideal interface, such as 
from the supposedly negatively charged LAO top sur- 
face (driven by the "polar catastrophe" ) fiS^i^ from oxy- 
gen vacancies f^ii^^—ii^ or from Sr-La intermixing at the 
interface . 

Here we note that the source of the carriers should not 
be in question. As explained by Fig. [l] the polar discon- 
tinuity at the STO/LAO interface provides an intrinsic 
source of electrons with a density of 0.5 electrons (e^) 
per unit cell (3.3x10^'' cm~^). In an ionic picture (based 
on Sr2+, La^+, Ti^+, and 0^~), LaAlOs can be regarded 
as a stacking of (LaO)"*" and (A102)~ planes along the 
[001] direction; LaO planes donate electrons which be- 
come bound in the neighboring AIO2 planes. SrTiOa, 
on the other hand, is composed of alternating charge- 
neutral (SrO)° and (Ti02)° planes. At the interface with 
Ti02-terminated STO, the Ti02 layer is already charge- 
neutral, causing the interfacial LaO layer to act as a sheet 
of donors, donating 0.5e^ per unit cell, which due to the 
large conduction-band offset between LAO and STO flow 
into the STO. The question is therefore not "where do the 
carriers come from?" , but rather, "where do the electrons 
disappear to?" I.e., why is the observed carrier density 
an order of magnitude lowe r"*^'^'^^!^^ than the expected 
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0.5e- per unit cell (3.3x10" cm-^)? 

Based on first-principles calculations combined with 
Schrodinger-Poisson (SP) simulations we attribute the 
problem to the lack of a suitable termination for the 
top surface of the LAO layer. The asymmetric nature 
of the resulting layer structure drains electrons away 
from the 2DEG, and further exposes the heterostruc- 
turc to the detrimental effects of point-defect forma- 
tion. Our insights enable us to propose specific strategies 
for overcoming the problems associated with the current 
STO/LAO structures, and also provide guidelines for the 
choice of other materials combinations. We note that 
2DEG densities as high as 3.3x10^'* cm~^ have already 
been observed at SrTiOs/GdTiOa interfaces, consistent 
with the predictions of our model. 

Our first-principles calculations are based on general- 
ized Kohn-Sham theory within the projector-augmented 
wave method as implemented in the VASP codei^"— We 
use the hybrid functional of Heyd, Scuseria, and Ernzer- 
hof (HSE)^ which has been shown to produce electronic 
structure and band gaps in closer agreement with exper- 
iment and a better description of charge localization,^^ 
compared to conventional functionals based on the local 
density approximation (LDA) or the generalized gradi- 
ent approximation (GGA). For the SrTiOa and LaAlOa 
cubic perovskite phases we use a 4x4x4 special k-point 
set for integrations over the Brillouin zone and an en- 
ergy cutoff of 500 eV for the plane- wave basis set. For 
the heterostructures, the calculations were performed for 
a (STO)8.5/(LAO)7.5 superlattice containing two equiv- 
alent Ti02-LaO interfaces, with a 4x4x1 k-point set 
and a 400 eV energy cutoff. Tests using 6x6x1 and 
8x8x1 meshes resulted in changes of less than 0.1 eV in 
Fermi energy, and charge densities remained essentially 
unchanged. The in-plane lattice constant was fixed to 
that of STO, representing a heterostructure coherently 
grown on an STO substrate. Full relaxation was allowed 
for both the out-of-plane lattice constant and all atomic 
positions. 

The SP simulations were performed using the 
nextnana^ simulation software, which solves for the elec- 
trostatic potential, charge density, and Fermi level across 
the heterostructures. The input parameters include elec- 
tron effective masses (1.0 me, fitted to reproduce the 
first-principles density of states) and dielectric constants 
(300 for STO, 27 for LAO). The valence-band band offset 
was set to the first-principles calculated value of 0.1 eV 
(higher in LAO). Background doping at a level n=10^^ 
cm~'^ in the STO was assumed. In the SP simulations we 
do not distinguish between the different types of carriers 
contributing to the 2DEG density, i.e., if d^y or dyz/dxz- 

The first-principles calculations quantitatively confirm 
that the free carriers in the 2DEG at the interface orig- 
inate from the interfacial LaO plane. Figure [5] shows 
the integrated charge density, plotted along the [001] di- 
rection, associated with the electrons generated at the 
interface. This charge density appears exclusively on the 
STO side and corresponds to occupied subbands in the 
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FIG. 2: (Color online) First-principles and Schrodinger- 
Poisson results for a (SrTi03)8.5/(LaA103)7.5 superlattice 
with Ti02-LaO interfaces, (a) Atomic structure, with oxygen 
atoms shown in red. Ti-centered octahedra are shown, (b) 
First-principles planar and macroscopically averaged charged 
density of the occupied subbands (excluding the STO and 
LAO valence bands), (c) Schrodinger-Poisson simulations for 
the same superlattice, showing good agreement with the first- 
principles results. CBM stands for conduction-band minimum 
and VBM for valence-band maximum. 



conduction band. The integrated charge is 3.3xl0~^^ 
cm~^ per interface (i.e., 0.5e~ per unit cell area), ex- 
actly what we expect based on the consideration of the 
interfacial LaO as a delta-doped donor layer. We observe 
that this charge is delocalized over multiple planes of Ti 
and not localized on a single Ti layer at the interface. If 
the latter were the case;^ the electrons would be immo- 
bilized on interfacial Ti^"*", which would be inconsistent 
with the observation of a high-mobility 2DEG. 

The symmetric nature of the charge density in the STO 
layer arises from the fact that our first-principles calcu- 
lation needs to maintain periodicity along the [001] di- 
rection, and therefore corresponds to a superlattice con- 
taining two interfaces. We have verified, however, that 
the results described here do not depend on the thickness 
chosen for the STO and LAO layers. 

Very similar results are produced by SP simulations^! 
for the (STO)8.5/LA07.5 superlattice. The nominal 
charge density of 0.5e^ per unit cell area, as obtained 
from first-principles calculations, is reproduced in the 
SP simulation for the (STO)8.5/LA07.5 superlattice, as 
shown in Fig. [D^c). While these simulations do not in- 
clude the intricacies of the STO conduction-band struc- 
ture, they do accurately capture the overall carrier den- 
sity and distribution near the interface, which is the fo- 
cus of our study. The SP simulations allow modeling 
systems with larger dimensions and, more importantly. 
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FIG. 3: (Color online) Schrodinger-Poisson simulations for 
SrTiOa/LaAlOs interfaces. Layer structures are depicted on 
the left, and the corresponding band diagrams and charge 
density distributions on the right. The zero of energy is placed 
at the Fermi level, (a) STO/LAO/STO with two equiva- 
lent Ti02-LaO interfaces. The integrated electron density 
is 3.3x10" cm"^ per interface, (b) STO/LAO/STO with in- 
equivalent interfaces: Ti02-LaO on the left, A102-SrO on the 
right, (c), STO/LAO with Ti02-LaO at the interface and 
an A102-terminated surface containing acceptor-like surface 
states (green horizontal bar), (d) STO/LAO with Ti02-LaO 
at the interface and a passivated surface. 



lacking periodicity, and hence enable us to study layer 
structures that are beyond the capabilities of the first- 
principles calculations. 

Given that our simulations so far indicate that the 
2DEG density should be expected to correspond to 0.5e~ 
per unit cell, we now address the question of why experi- 
ment shows much lower densities.— >^'i^'i^ The answer lies 
in the fact that the type of symmetric structure depicted 
in Fig. [21 with two identical interfaces, is never achieved 
experimentally. In most practical implementations, the 
LAO layer is of finite thickness and has a surface termi- 
nated on an AlO 2 plane. The consequences are examined 
in Fig. |3] for a 2-nm thick LAO layer (about 5 unit cells) 
on thick STO. In panel a it is assumed that a perfect STO 
layer (with Ti02-LaO planes at the interface) can be de- 
posited on top of the LAO. This effectively reproduces 
the symmetric situation that was investigated in Fig. [2l 
the only difference being the smaller LAO thickness and 
larger STO thickness; reassuringly, the results are very 
similar, with the full density corresponding to 0.5e~ per 
unit cell appearing in the 2DEG. Note that the electro- 
static potential (reflected in the slope of the conduction 
band) is essentially flat across the LAO layer, indicating 
no charge is being transferred between the two interfaces. 



The case depicted in Fig. [3l^a) is unrealistic, because 
experimentally it has turned out to be difficult (or even 
impossible) to grow STO/LAO/STO structures with an 
interface between LaO and Ti02 planes on the right- 
hand sideJ^ Figure [Sljb) depicts the situation for an 
STO/LAO/STO layer structure with an interface be- 
tween AIO2 and SrO on the right. By similar logic as 
applied to the Ti02-LaO interface leading to donor dop- 
ing, an A102-SrO interface leads to acceptor doping with 
a sheet density of 3.3x10"^'' cm~^ since the AIO2 layer 
is lacking 0.5e~ per unit cell which (in the bulk) would 
come from an LaO plane (the SrO plane in STO being 
charge- neutral) . The electrons at the Ti02-LaO inter- 
face are higher in energy than the holes at the A102-SrO 
interface, making this situation unstable. The energy of 
the system is lowered by transferring electrons from the 
Ti02-LaO interface to the A102-SrO interface, leading to 
a strong dipole across the LAO layer (note the slope in 
electrostatic potential). In the process, the 2DEG den- 
sity is drastically reduced, although not to zero because 
the electric field in the LAO layer results in the valence 
band of LAO rising above the Fermi level (with a hole 
gas appearing at the right-hand interface), thus limiting 
the transfer of electrons. 

Having an ideal A102-SrO interface on the right 
[Fig. [3Jb)] is similar to terminating the LAO layer with 
an ideal AIO2 surface — indeed, all planes in the STO on 
top are charge-neutral and hence do not contribute to any 
charge exchange. The only difference is that a realistic 
A102-terminated LAO surface would exhibit oxygen dan- 
gling bonds, giving rise to partially filled surface states 
with energies in the lower part of the band gap. Simi- 
lar charge transfer would occur as described for the case 
depicted in Fig. [3]^b) , except that the acceptor states are 
now deep acceptor states (with an ionization energy of 1 
eV above the LAO VBM), as shown in Fig. [S^c). The 
Fermi level at the surface will be pinned at these accep- 
tor states, limiting the rise of the LAO VBM. This is 
actually beneficial since it limits the amount of electron 
transfer out of the 2DEG, compared to the situation of 
Fig. El^b). Note that this corresponds to a decrease in 
the slope of the potential, i.e., the magnitude of the elec- 
tric field, across the LAO layer: the smaller the field, the 
higher the 2DEG density. 

This suggests a strategy for increasing the electron 
density in the 2DEG, namely minimizing the slope in 
the potential across the LAO. Figure [3l^d) depicts a fully 
passivated surface, i.e., the density of acceptors at the 
surface is assumed to be zero. While there is a still a 
slight slope in the potential (due to the fact that the cen- 
ters of gravity of the positive and negative charge distri- 
butions do not coincide), the 2DEG density now recovers 
its nominal value of 0.5e~ per unit cell. 

Our findings easily explain the experimentally ob- 
served dependence of 2DEG density on LAO thickness.^- 
As noted in the discussion of Fig.[3^b) above, full transfer 
of electrons to the top surface will tend to occur if the 
LAO layer thickness is insufficient to bring the VBM at 
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the surface above the Fermi level. Given a positive sheet 
charge density of 3.3x10"" cm'^ at the STO/LAO in- 
terface, Gauss' law predicts a field of about 0.25 V/A, 
and thus it takes a "critical thickness" of about 3 or 4 
unit cells of LAO to develop enough of an increase in po- 
tential to bring the VBM (or acceptor-like surface states) 
of LAO above the Fermi level, at which point transfer of 
electrons out of the 2DEG is suppressed and observable 
mobile charge appears in the 2DEG. This is confirmed 
by explicit SP simulations as a function of LAO layer 
thickness (not shown). In principle the carrier density 
in the 2DEG should continue to rise as a function of 
LAO thickness. However, growing thicker LAO layers 
of high quality will be difficult because an increasingly 
large fraction of the layer has the Fermi level lying close 
to the VBM, a condition that is conducive to formation 
of oxygen vacancies. These will both interfere with the 
charge balance and degrade the quality of the layer. All 
of this is completely consistent with results in the lit- 
erature; our main point is that complicated arguments 
(e.g., relating to interfacial reconstructions) are unnec- 
essary to explain the experimental observations If the 
exposed AIO2 surface is passivated, our model predicts 
that electron transfer is suppressed and a 2DEG can in 
principle be observed, even for LAO layers below the crit- 
ical thickness. This explains why some experiments have 
observed conducting interfaces for LAO thickness of less 
than 4 unit cells 

We now turn to the observed variation in 2DEG 
density with oxygen partial pressure during growth or 
annealing ji^ii^ LAO layers are typically A102-terminated 
and exhibit partially filled oxygen dangling bonds. Given 
the high electronegativity of oxygen and the position of 
such dangling-bond states close to the VBM of LAO, a 
strong driving force exists to fill these surface states with 
electrons. Modifications of the surface that remove oxy- 
gen dangling bonds suppress electron transfer to the sur- 
face, and hence lead to a higher 2DEG density. Growth 
or annealing in an environment with low O2 partial pres- 
sure results in surface reconstructions containing oxygen 
vacancies (or equivalently cation adatoms) , effectively re- 
moving oxygen dangling-bond states from the surface. 
The situation is then closer to the scenario of the pas- 
sivated surface described in Fig. ^d), which shows that 
a high-density 2DEG can develop. While such oxygen 
treatments cannot be expected to lead to full passiva- 
tion, the predicted trend of lower oxygen pressure re- 
sulting in higher 2DEG density is definitely consistent 
with experiment, without having to invoke modification 
of or point-defect formation at the buried STO/LAO in- 
terface. Indeed, we consider the latter unlikely due to 
the high formation energy of oxygen vacancies in n-type 
STOJi 

Conversely, annealing under high O2 partial pressure 
leads to a higher density of oxygen-related surface states 
that will consume electrons from the 2DEG. This ex- 
plains the seemingly puzzling fact that attempts to per- 
fect the structural quality of STO/LAO heterostructure 



by growing or annealing under high O2 partial pressure 
often lead to high sheet resistance or insulating behavior 
at the interfacei^'S'i^ 

Manipulating oxygen partial pressure may not be the 
most effective means of passivating the LAO surface. 
Hydrogen tends to be a good passivating agent. In- 
deed, first-principles calculationsSS found that hydro- 
genation of the LAO surface leads to an increase in the 
2DEG density, although we disagree with Son et alJs 
interpretation^ that hydrogen donates electrons to the 
interface. The correct picture, as argued above, is that 
hydrogen passivates the surface and in the process pre- 
vents electrons being drained away from the interface. 
Experimental efforts to identify the most effective means 
of controlling and passivating LAO surfaces could be 
highly fruitful. 

Another approach to prevent electron transfer to the 
surface is to provide a source of electrons to the surface, 
for instance by depositing a metal on top of the LAO. 
The effect will depend on the metal used, specifically, 
on the work function of the metal relative to the elec- 
tron affinity of STO. Metals with work function larger 
than the electron affinity of STO (e.g., Au) will still re- 
sult in suppressed 2DEG densities since they will not 
succeed in suppressing electrons draining away from the 
interface. Metals with work functions equal to or smaller 
than the electron affinity of STO (such as Ti or Al) are 
needed to increase 2DEG densities. These qualitative in- 
sights are confirmed by explicit SP simulations. These 
effects have also been observed in recent first-principles 
calculations}^ although again the interpretation should 
not be that electrons are being transferred from the metal 
to the interface. Ultimately, though, metal capping layers 
on STO /LAO heterostructures may be more of academic 
interest, since it may prevent experimental probing of 
the 2DEG and also prove incompatible with device ap- 
plications. Emphasis on surface passivation techniques, 
as described above, is a more promising route. 

Finally, we note that our model of the fundamental 
physics at STO/LAO interfaces transcends the specific 
materials system being discussed here and is generally 
applicable to oxide interface. For instance, it explains 
why a 2DEG with the full nominal density of 3.3xl0~^* 
cm~^ has been observed at SrTiOs/GdTiOs interfaceSf2S 
GdTiOa (GTO) is composed of alternating positively 
charged (GdO)"*" planes and negatively charged TiO^ 
planes, similar to LaAlOs (note that Ti has valence 3 
in GTO). An STO/GTO interface with Ti02-GdO inter- 
facial planes will therefore also act as a sheet of donors. 
The difference with the STO/LAO case lies in the fact 
that no electric field occurs within the GTO layers. In- 
deed, STO can be grown with high quality on top of 
GTO, and since the interfaces are always between Ti02 
and GdO planes the top and bottom interfaces of each 
GTO layer are identical by construction. This symme- 
try prevents an electrostatic potential buildup, as shown 
in Figs. H and [3)^a) and allows the fufi 3.3x10"" cm"^ 
density to be present in the 2DEG. 



5 



Even in the absence of an STO overlayer, a buildup 
of electrostatic potential is unlikely in the GTO layer. 
If GdO-terminated, both the interface and surface would 
exhibit donor-like behavior. If Ti02-terminated, the elec- 
trons that are needed to fill acceptor states do not flow 
into deep-lying oxygen-derived states, but rather into a 
Ti-derived lower Hubbard band which lies not far below 
the CBM of This position of the Fermi level at 

the surface again suppresses a buildup of potential and 
maintains the full 2DEG density at the STO/GTO inter- 
face. 

In summary, based on first-principles calculations and 
Schrodinger-Poisson simulations we demonstrate that 
electronic conductivity at the STO/LAO interface arises 
from electrons that are intrinsic to the interface. This 
precludes the need to invoke other sources of electrons 
such as the top LAO surface (according to the po- 
lar catastrophe model) ji^i^ or oxygen vacancies acting 
as donors.- The suppression of the 2DEG density at 
STO/LAO interfaces has often been attributed to inter- 
facial reconstructions (either atomic or purely electronic, 
based on mixed valence of Ti), which in turn were in- 



voked as a consequence of a "polar catastrophe". Our 
present results show there is no need for invoking such 
mechanisms. They emphasize the need for measures to 
prevent electrons draining away form the interface, which 
can be accomplished by preventing an electrostatic po- 
tential buildup in the LAO layer. Proposed strategies 
include passivation of the surface, or depositing metals 
with suitably low work functions. These insights into the 
origin of carriers at the STO /LAO interface will pave the 
way to enhanced control of the 2DEG at the interface of 
complex oxides. 
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